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Fermionic Molecular Dynamics (FMD) is a microscopic many-body approach which has been 
used successfully to study the structure of light nuclei in the p- and ^J-shell. FMD uses a Gaussian 
Q ■ wave packet basis which contains the harmonic oscillator shell model and Brink-type cluster 

^ . model wave functions as limiting cases. A realistic effective interaction that has been derived from 

the Argonne VI 8 interaction by treating explicitly short-range central and tensor correlations is 
employed. 

We present here a first application of the FMD approach to low-energy nuclear reactions, namely 
' the ^He(a,7)^Be radiative capture reaction. We divide the Hilbert space into an external region 

■ where the system is described as ^He and ^He clusters interacting only via the Coulomb inter- 

action and an internal region where the nuclear interaction will polarize the clusters. Polarized 
configurations are obtained by a variation after parity and angular momentum projection proce- 
dure with respect to the parameters of all single particle states. A constraint on the radius of 
the intrinsic many-body state is employed to obtain polarized clusters at desired distances. The 
boundary conditions for bound and scattering states are implemented using the Bloch operator. 
The FMD calculations reproduce the correct energy for the centroid of the 3/2~ and l/2~ bound 
. states in ^Be. The charge radius of the ground state is in good agreement with recent experimental 

results. The FMD calculations also describe well the experimental phase shift data in the 1/2+, 
3/2+ and 5/2+ channels that are important for the capture reaction at low energies. Using the 
bound and scattering many-body wave functions we calculate the radiative capture cross section. 
The calculated S factor agrees very well, both in absolute normalization and energy dependence, 
with the recent experimental data from the Weizmann, LUNA, Seattle and ERNA experiments. 
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1. Introduction 



Low-energy nuclear reactions play an important role in many astrophysical scenarios. In many 
cases experimental data are not available at the energies relevant for the astrophysical processes. 
Typically such reactions are described with potential models which describe the reaction partners 
as point-like nuclei interacting via nucleus-nucleus potentials fitted to experimental data on bound 
and scattering states. In a microscopic ab initio picture however, the system is described as a many- 
body system of interacting nucleons. The wave functions are fully antisymmetrized and realistic 
nucleon-nucleon interactions are used. 

In the Fermionic Molecular Dynamics (FMD) approach we aim at a consistent description 
of bound states, resonances and scattering states using realistic low-momentum nucleon-nucleon 
interactions obtained in the Unitary Correlation Operator Method (UCOM) [^. Intrinsic many- 
body basis states are Slater determinants using Gaussian wave packets as single-particle states. This 
basis contains harmonic oscillator shell model and Brink-type cluster wave functions as special 
cases. The symmetries of the system are restored by projection on parity, angular momentum 
and total linear momentum. The many-body eigenstates of the realistic Hamiltonian are obtained 
in multiconfiguration mixing calculations. FMD has already been used successfully to describe 
nuclei in the p- and ^J-shell, like the Hoyle state in ^^C [^] or halo- and cluster-structures in Neon 
isotopes [0]. 

In this contribution we present results for the ^He(a,7)^Be capture cross section using the 
FMD approach. This reaction has been studied using potential models [§, ^] and microscopic clus- 
ter models [^, |8[] assuming ^He+'^He cluster wave functions. Polarization effects were considered 
by including the ^Li+p channel in [Q, [l^. Consistent ab initio calculations starting from realistic 
interactions have not been possible up to now. First attempts using Variational Monte-Carlo [ pi] ] 
and the No-Core Shell Model [|12|] only calculated the asymptotic normalization coefficients from 
the bound state wave functions. 




Figure 1: Density distributions of the intrinsic basis states. Top: frozen cluster configurations, bottom: 
selected polarized configurations obtained by variation after angular momentum projection. 
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2. Fermionic Molecular Dynamics 



In FMD the intrinsic many-body basis states are Slater determinants 

|e> = ^{|^i>«)...«) |^A>}, 

using Gaussian wave packets as single-particle states 



(2.1) 



(2.2) 



The complex parameters bk encode the mean positions and momenta of the wave packets, the width 
parameters Uk are variational, the spins can assume any direction. To restore the symmetries of the 
Hamiltonian the intrinsic wave function is projected on parity, angular momentum and total linear 
momentum P = 0. 

|(2;/^M7r;P = 0> =Pi^^P^P^=^|e>. (2.3) 

As in a microscopic cluster model we have "frozen" cluster configurations where "^He and 
^He clusters are located at a distance R. These Slater determinants can be rewritten using RGM 



basis states [|13|, [14j]. The RGM representation can be used to match the logarithmic derivative of 
the relative wave function of the clusters to the asymptotic Whittaker or Coulomb wave functions. 
In the interaction region we have additional FMD basis states that are necessary to describe the 
polarization of the clusters. These basis states are generated by variation of all single-particle 
parameters after parity and angular momentum projection. We use a constraint on the radius of the 
intrinsic states to generate configurations corresponding to different cluster distances. In Fig. [l| we 
show intrinsic density distributions for some frozen and polarized basis states. 



3. Bound and Scattering States 



Using the microscopic /^-matrix approach of the Brussels group [|13|, |14j] we solve the Schrodinger 
equation with the proper boundary conditions for bound and scattering states. With frozen configu- 
rations alone the 3/2~ and l/2~ states are bound by only 240 keV and 10 keV, respectively. In the 




E [MeV] 




E [MeV] 



Figure 2: ^He-^He scattering phase shifts for P-waves left and the S and Z)-waves right. Dashed lines show 
results using only frozen cluster configurations, solid lines show results for the full model space. 



3 



Microscopic Nuclear Structure and Reaction Calculations in the FMD Approach 



Thomas Neff 






Hi88 
Kr82 




Os82 




Pa63 


• 


NS04 


• 


Co07a 


• 


Co07p 


• 


Br07a 


• 


Br07p 


• 


Le09r 


• 


Le09a 




De04 




FMD 



1 1.5 
Ecm [MeV] 



Figure 3: The astrophysical ^'-factor for the ^He(a,7)^Be reaction. The FMD result is given by a black solid 
line. Old data are shown as gray symbols together with an 7?-matrix fit. Recent data are shown by colored 
symbols. 



full Hilbert space, including the polarized configurations, we obtain binding energies of 1.49 MeV 
and 1.31 MeV with respect to the cluster threshold. We therefore reproduce the centroid energy 
but underestimate the splitting of the 3/2~ and l/2~ states. This will effect the branching ratio but 
not the total cross section which essentially only depends on the centroid energy. The calculated 
charge radius of 2.67 fm is in good agreement with the experimental value of 2.647(17) fm [p^]. 
In Fig. ^ we present the calculated phase shifts together with the experimental data [|T6[]. Again 
we see a sizable effect when we compare the results using only the frozen configurations with the 
results obtained in the full Hilbert space including polarization effects. 



4. Capture Cross Section 

Using the microscopic bound and scattering wave functions we calculate the radiative capture 
cross section. We have contributions from the internal region calculated with the FMD many- 
body wave functions and from the external region up to large cluster distances calculated with the 
matched Whittaker and Coulomb solutions. At low energies we can restrict ourselves to electric 
dipole transitions from the S- and D-wave channels. The obtained cross section, presented in form 
of the astrophysical S factor is shown in Fig. ||. It agrees very well both in absolute normalization 
and energy dependence with the recent experimental results obtained at the Weizmann Institute 
[Q, the LUNA collaboration [|l|], at Seattle [jl^ and by the ERNA collaboration [0]. 
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